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ABSTRACT 

The magnetorotational instability(MRI) in the Kerr spacetime is studied on 
a 3+1 viewpoint. The Maxwell's equations are expressed in a circularly orbiting 
observer's frame that co-rotates with matter in Keplerian orbits. The hydromag- 
netic equations are represented in a locally nonrotating frame(LNRF). There 
exist large proper growth rates in MRI around a rapidly rotating black hole. The 
large " centrifugal force" and the rapid variations of magnetic fields are caused by 
the rotation of spacetime geometry. As the result, in the extreme Kerr case the 
maximum proper growth rate at r = r ms becomes about twelve times as large as 
that in Schwartzshield case, where r ms is the radius of a marginally stable orbit. 
The unstable range of wavenumber expands according to the rotational speed of 
spacetime geometry. The over stable mode of the instability becomes remark- 
able when the circular motion of the disk is quasi-relativistic and the strength of 
magnetic field is so large that the Alfven velocity is va > 0.1c. When the waves 
of perturbations propagate in the radial direction, these waves oscillate and their 
amplitudes grow exponentially. This instability is caused by the differential rota- 
tion of the circularly orbiting frame. The universality of the local Oort A-value 
of the disk is discussed in curved spacetime. In the extreme Kerr geometry, the 
amplitude of the maximum growth rate in the dynamical shear instability (DSI) 
reaches to infinity at r = r ms . The behaviors of maximum growth rates between 
MRI and DSI are remarkably different each other. 

Subject headings: accretion — hydromagnetics: instabilities — turbulence: rela- 
tivity — black hole 



1. Introduction 

The various intensive phenomena, x-ray or 7-ray emission, jet, etc. are found in many 
active objects. Some intense 7-ray bursts (GRBs) are believed to arise from the hot, dense 
matter accretion onto a rotating black hole (Popham, Woosley Sz Fryer 1999; Di Matteo, 
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Perna, & Narayan 2002; Kohri, & Mineshige 2002; Yokosawa, Uematsu & Abe 2004). The 
viscosity is the important physical quantity to balance the heat in the disk and to determine 
the dynamics of the accretion. The magnetohydrodynamical turbulent viscosity is driven by 
the magnetorotational instability (MRI) (Balbus & Hawley 1991; Hawley & Balbus 1991; 
Hawley, Gammie, & Balbus 1995; Brandenburg et al. 1995; Matsumoto & Tajima 1995; 
Stone et al. 1996; Balbus & Hawley 1998). We study here the MRI on the basis of the 
general relativistic magnetohydrodynamics. 

The MRI in the Kerr metric had been considered by Gammie (2004) . He had analyzed 
the dynamical shear instability(DSI) for circularly orbiting particles. The maximum growth 
rate in MRI was then discussed according to the analogy between DSI and MRI. It had been 
previously shown by Balbus & Hawley (1992) that in a Newtonian case the both maximum 
growth rates in MRI and DSI coincide each other. However its consistency is not always 
satisfied in general spacetime. The MRI in curved spacetime should be investigated on the 
general relativistic Maxwell's equations and hydrodynamics, and then be examined for the 
consistency of the maximum growth rate between the two types of instabilities. 

The generation of MRI is caused by the two processes; the strengthen of magnetic 
field due to the differential rotations of the disk, and the transfer of angular momentum by 
magnetic stress (Balbus & Hawley 1998). The rotation of the spacetime geometry exerts the 
characteristic effects on the kinematics and on the induction of magnetic field. In the rotating 
spacetime geometry, there appears "the gravitomagnetic force" which is an inertial force in 
3+1 point of view (Thorne et al. 1986; Punsly 2001). The gravitomagnetic force depends 
on the gradient of the geometrical angular velocity, Vw. The gravitomagnetic force exerts a 
torque on a spin outside a black hole (Thorne et al. 1986) and strengthens the centrifugal 
force. The increased radial velocity then induces the large variation of magnetic fields which 
transfer the angular momentum. The rotating velocity of the disk becomes extremely large 
in the vicinity of a rapidly rotating black hole. The changing rate of magnetic field due to 
the differential rotation of the disk is then largely increased. The generation of MRI may be 
very active around a rapidly rotating black hole. 

It is useful to introduce the 3+1 formulation of general relativity which is excellent to 
make clear the physics on the stability. The evolutions of the perturbations are then simply 
evaluated by the analogous terms to the Newtonian dynamics and to the Maxwell equations 
in flat spacetime. We introduce the Keplerian orbiting observer's frame (orbiting frame) to 
evaluate the electromagnetic fields perturbed by the motion of fluid. On the other hand 
the dynamic equations of fluid interacting with magnetic fields are analyzed in the locally 
non-rotating frame (LNRF)(Bardeen et al. 1972) since this frame becomes an inertial frame 
at the far distance from a black hole. Thus the dispersion relation expressed in these frames 
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approaches to a Newtonian one at the far distance. 

When the Maxwell equations are described in the orbiting frame, there appear explicitly 
the electromagnetic fields generated by the frame moving . When the rotating velocity of 
the frame is relativistic, the induction of magnetic field by the differential rotation of the 
orbiting frame becomes efficient. This generation of magnetic field causes an overstable mode 
of instability. We investigate here the characteristic behaviors of this new mode in MRI. 

The intensive phenomena in active objects may be produced by the nonstationary mag- 
netohydrodynamical processes in the accretion disk around a rotating black hole. We are 
interesting in the spatial distribution of a growth rate of MRI. Yokosawa, Uematsu & Abe 
(2004) had represented the angular velocity of circular orbits for a geometrically thick disk, 
fl(r,z). By using this angular velocity Q(r,z), we investigate the distribution of the maxi- 
mum growth rate T max (r, z)over the space around a rotating black hole. 

When the geometry of spacetime rotates, it is important to select the suitable frame for 
the analysis of the instability. In the Boyer-Lindquist coordinate frame (BLCF) the dragging 
of inertial frames becomes so severe that the t coordinate basis vector (d/dt) goes spacelike 
at the static limit tq (outer boundary of the ergosphere). In LNRF the frame dragging effects 
of the hole's rotation are canceled out. Gammie (2004) investigated the DSI in BLCF. 
We reanalyze the upper bound to the growth rate of DSI in LNRF by using the general 
relativistic Hill equations. We then examine the asymptotic behaviors of the growth rates 
for r — > r ms when the frame is exchanged from BLCF to LNRF. We inquires the suitable 
frame to evaluate the instability, MRI or DSI. The universality of local Oort A-value in the 
shear instability is then examined by comparing the dispersion relations for DSI and MRI 
in the suitable frame. 

In the next section we shall describe the Maxwell's equations in the orbiting frame and 
the magnetohydrodynamic equations in LNRF. The dispersion relation and the properties 
of growth rates in curved spacetime are analyzed in the section 3. The universality of the 
growth rate of the shear instability in general, gravitational field is discussed in section 4. 

2. Basic Equations 

The preferred frame describing the basic equations is selected definitely to evaluate the 
perturbations. The LNRF becomes an inertial frame at the infinite distance from a black 
hole. Thus the characteristic frequencies in dynamics, e.g., epicyclic frequency, are definitely 
expressed in LNRF. The perturbations of hydrodynamics are represented in LNRF. The 
perturbations of magnetic field exerting on the fluid are simply evaluated in the frame co- 
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moving with the fluid in a Keplerian orbit (Punsly 2001). The both frames, LNRF and 
orbiting frame, are combined with the Lorenz transformation. 

We choose units with G = c = 1. The metric of spacetime, in terms of Boyer-Lindquist 
time t and any arbitrary spatial coordinates x\ has the form 

ds 2 = -a 2 dt 2 + g jk {dx J + (3 j dt){dx k + fi k dt). (1) 

The metric coefficients a, ft, are the lapse, shift, and 3-metric functions. Theses functions 
are given by 

/£A\i/2 E sin 2 9 A 

« - y—^- J , 9rr - 9ee - S, g w - — - — 

2 Mar 

fif = —uj = — — , ft = gj k = for all other j and k, 

where M is the mass of the black hole, a is its angular momentum per unit mass, and the 
functions A, E, A are defined by 

A = r 2 -2Mr + a 2 , E = r 2 + a 2 cos 2 9, A = (r 2 + a 2 ) 2 - a 2 A sin 2 9. 

We introduce a set of local observers who rotate with the Keplerian angular velocity. 
Each observer carries an orthonormal tetrad. For the Keplerian orbiting observer with a 
coordinate angular velocity of Q, its world line is r =constant, 9 =constant, ip = flt+ 
constant. The observer in the LNRF who rotates with the angular velocity u measures the 
linear velocity of a particle moving in a Keplerian orbit as = dx^/dx^ = a~ 1 ^/g lfilfi (Vl — 
uj) where a;*- -* and x^ are the time and spatial coordinates in LNRF. Its Lorentz factor is 
7 = (1 — v 2 ^)^ 1 / 2 . We adopt the orbiting frame with the set of it's basis vectors: 



_ d _ 2 d 

where e", e^, e M1 , e' 12 are difined by 

e» = a , e 2 ^ = g w , e 2 ^ = g rr , e 2 ^ = g ee . (4) 

The corresponding orthogonal basis of one-forms is 

uj 1 = adt - v^e^dip, uj* = -Q^dt + -fe^dip 

uf = e w dr, uf = e^d9. (5) 
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Here a is defined by a = 0:7(1 + v^v^), where v w is the "rotating velocity of the geometry 
of spacetime", v w = uoe^~ u . The rotation one-forms, a;?, which allow one to read off the 
connection coefficients by a;? = ^yp l . Hereafter the components of the tensor in the 
orbiting frame are depicted by the superscript with hat. The connection coefficients r~ are 
represented in Appendix A. 

The Maxwell equations for the electromagnetic field F^ v produced by a four current 
density are 

/•':'/: = -I"- Fptt + 1\,,„ + = 0- (6) 

In the orbiting frame they are represented as follows: 

^ + -Vx(oE) 

or a 

- n^e%B- VlnQ-e b v {v>) • BV In d - x Vln(fi/d) = 0, (7) 



^ - -Vx(oB) 

or a 



- n 



^e%E • V In n - e b v {tp) • EV In a - 7 \ v) vqB x V In (SI/ a) = — —J, (8) 

Air 

V-E - t; M E-Vn + e\ ) -(BxV)lnd = -J , (9) 
V B - v {v) B • Vfi - e b v M ■ (E x V) In a = 0. (10) 

Here rff is the variation of time in the orbiting frame and it combins with the coresponding 
time dt in the static frame by df = j^adt. The symbol is the unit vector in the % direction 
in the space of the orbiting frame. The vectors B = B % w.i and E = E 1 ^ are the magnetic 
and electric fields with components E r = F 0r , = F re , and so on. The symbol V is the 
covariant derivative operator, and vy^,), vq and e b are defined by V( v ) = iv^n^, vq = i\e^~ v 
and e b = 7 2 (l+t> (u,)vj). When the angular velocity of the frame is taken to be the geometrical 
angular velocity, Vt — > u, the above expressions of Maxwell equations (7) — (10) result in the 
forms described in LNRF (Thorne & Macdonald 1982; Thorne et al. 1986). 

The Lorentz force is represented by 

1 



J x B = -V 

c 



f§- + i B ' V(dB) + f W „(-|-V + |B ■ V) ln(d/!2). (11) 



The energy-momentum tensor in LNRF, T^^\ is expressed with the four velocity of 
fluid the pressure P, and the density p as 

T m) = ( p + p) u « u (j) + rj(m P _ ( 12 ) 
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Hereafter the components of a tensor in LNRF are expressed by the superscript with the 
parenthesis. The tensor 77WO') is the Minkowskian metric. The Euler equations with the 
Lorentz force, F^^J^, are represented by 

d(p + + ^ a{p + p)uu(i) + l^ ap + pB^M = F (0M j (i/)> (13) 

where f and V are the time and the covariant derivative operator in LNRF. The variation 
of time, df, is combined with df by df = jdf. The Lorentz force is related likewise as 
FWM J (l/) = 7F Jl> J,>. The term rg )(l/) TWM is the "inertial force" in the (i) direction. The 
forces for (i) = (r), (0) are expressed in explicit forms: 

r S)M TMM = ^ ( ° )(0) V (l) In a + T(°)^u w V« In w - T^^ V (i ty 

+ T«^V°V - T^V^nj. (14) 

Here the (r) component of T^^V^Vi -T^^VW/ij is T^W V^i -T^^VW^, and the 
(6>) component is T^W V (r) /i2 - T (r)(r) V w /xi- Th e (</?) component of the "inertial force" is 

rg (l/) TMM = (p + p)«^u.V^. (15) 

This term is correspond to the "Coliori's force". 
The equation of continuity is 

dpu^ 1~ 

+ -V • apu = 0, 16 

or a 



3. Groth Rates of Magnetorotational Instability 

3.1. Dispersion Relation 

The magnetohydrodynamical variables are expanded by the infinitesimal quantities, 
v = v + 5v, B = B + <5B, p = p + 5p and P = P + SP. We assume the axisymmetric 
perturbations and use the local WKB approximation. The variables are unified to be the 
components in the orbiting frame by the transformations of the hydrodynamic equations 
from LNRF to the orbiting frame. The linearized equations are then described by 

V-5v = 0, (17) 
d5v f ~ r & 1 85 P bp dP 

— - - 2Qf6v v + — — — ~r ~z ~ 

or 7p 8x r 7/tr 8x r 
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+ 



B§ d6B<> B§ d5B* 



4irp V dx 



dx 6 



+ 2B°5B° 1 \ tp) v a V r In (a/ft) = 0, 



d5v e 

~~dT 

dSv^ 
~dr~ 
d5B' 

dr 
d5B § 

~dr~ 

d5B* 

~dr~ 

5 dSp 
3p df 



+ K ■ 5v 



Airp 
B • V5v f = 0, 



IP dx IP 1 dx 9 
1 = 0, 



B • V5v e - B § 5v f 7 2 v {tp) v n V f In (a/ft) = 0, 
B • - 7 2 wn(5B • V In ft = 0, 



+ 8v* 



din Pp~ 5 / 3 e ndlnPp 
— -. h for 



-5/3 



0, 



where the coefficients ft = (ftf , ft#) and k = (k?, Kg) are defined by 



ft = 7(%>)^ — TT^Vlncu), and 



K 



= 7^)Vln (7 2 e^ (¥)) ). 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 

(25) 



With the spacetime dependence of perturbations, e % ( k t xT + k e x " Tt ) ; the explicit dispersion 
relation is derived as 



B 



K + ( 7^ - A^ r 



2n 



r 2 -(« 2 + x 2 )(k.v^ 



+ ^7 v iv) vn V r In (ft/a) 2 7 2 ^)^ V r In (ft/a) - i3fc* 



(k ■ v A ) 



-f 2 = 0, (26) 



where 



r^r^-(k.v^) 2 , e = kf + kl n 



3 <9P 9 In Pp~ 5 / 3 
57p <9x* 



for i = r,9, (27) 



and 



ft 



off)- - — O-Vlf- - — -- 2 - °- 2 - ^ kf -^\(^r kr 

\ r k )\ k~ /' 



X 2 = -2 7 2 ^ ( ftf - (V r - ^ V e J In ft. (28) 



3.2. Growth Rates 

The dispersion relation (26) has the resemblant form to a Newtonian case except the 
last term. When v \ <C N 2 = N~ + iV|, i.e., the magnetic pressure is negligibly small in 
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comparison with the gas pressure, j3 = P gaS s/ 'Pmag 1, the last term containing v\ r y 2 V( ip - ) vn 
can be ignored in (26). At the equatorial plane and far from a black hole, the functions, k 2 
and x 2 , have the following forms: 

r or o In r 

The equation (26) then coincides with a Newtonian dispersion relation (Balbus & Hawley 
1991). The equation (26) is the general dispersion relation of MRI for all radii and for all a. 

In the Newtonian case, the maximum growth rate relative to a angular velocity of a 
circular orbit, T max /Q, has an universal value in all radii and in all mass scales of a central 
object (Balbus & Hawley 1991). Here we introduce the angular velocity of a circular orbit 
measured in LNRF, = (Q — uo)/a. At the next it is transferred into the orbiting frame as 
ft = = (J} — u)j/a. It may be said that — Q is the proper angular velocity of the geometry 
measured in the orbiting frame. Let's consider first a simple case for /3 <C 1, k? = B r = 0, 
and introduce the dimensionless variables with use of a angular velocity & as 

f = T/n, k = (k • v A )/tl, N = Nf/Cl, Qf = ^fM R = K/Cl, and x = x/&- (30) 

Here it may be noticed that the dimensionless growth rate T coincides with that in LNRF,T = 
F/Cl = T/fl, where f is the growth rate measured in LNRF, f = —id\n5/df. The disk is 
assumed to be isothermal and the Brunt- Vaisala frequency is taken to be N~ = N?/Q 2 = 0.8 
and N? = N?/{l 2 = OMNf. 

The dispersion relation (26) will lead to instability for values of k less than the critical 
value: 

klit = f~ N 2 . (31) 

The unstable mode of a perturbation has a maximum growth rate, T maxi at the wavenumber, 
k — — kfyi&x such cts 

f - = -T (1 i+S >2 ' alld = f TtS§ (1 + (2R2 + R2)lf) - (32) 

When a = M, the asymptotic behaviors of the proper growth rate r ma rr(?") and wavenum- 
ber kmax{r) for r — > r ms are simply expressed by the asymptotic relation between Q(r) and 
oj{r). We show in Fig. 2 the variations of dimensionless variables, x 2 ^ 2 and ^ max ,k max , 
along the radius from a event horizon for a = M, 0. When r — > r ms and a — > M, the 
angular velocity of a circular orbit coincides with the geometrical angular velocity, Q — > uj. 
The "Coliori's force" n r 5v r in (20) then reduces to be zero, n r 5v r oc v ^ V r In (^y 2 e^v ^) = 
(Q — Ld)Ar~ 3 V r In (^y 2 e^v ^) — > for r — > r ms . Thus the epicyclic frequency, /t = (2fi f: ^) 1 / 2 , 
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becomes to be zero. Since x 2 > 2, the ratio N 2 /x 2 is negligibly small in (32). The maximum 
growth rate T max and wavenumber k max at r « r ms are then simply expressed as, 

f 2 ~ — — h 2 ~ ^— (11') 

1 max ~ ^ ft mai ~ ^ ■ l^J 

The value of x 2 is determined by the strengths of two actions. One is the radial accel- 
eration 5v r by the "centrifugal force", 5v r oc 2Qf5v^. Other is the induction of magnetic 
field by the differential rotation of fluid motion, 5B^ oc 8B r ^ 2 voV r In Q. We define the shear 
term by VL shear = r y 2 vaV r ln(fi). Then the function x 2 is expressed as x 2 — —^f^shear- 
The asymptotic forms of fl s hear and Qf are represented as fl s hear = — 1.5S1 and Cl? = for 
t — > oo, and Vt s h ear = — 4f2 and (if = —2~fV r uj = (A/\f?>)uj for r — > r ms . The efficiencies of 
the two actions relative to the angular velocity f2 or the proper angular velocity Cl(see Fig. 
2) are increased at r « r ms . When r — > r ms , the function fi^ in (25) is represented mainly by 
the variable related with the differential rotation of the geometry, v UJ V r lnuj, since another 
variable v ( v ,)V r ?/> reduces to be zero, V( lfi )'V r ip — (Jl — cj)e^~ iy V r -?/> = (fi — c<j)Ar _3 V r -?/' — > 0. 
The finite values of the growth rate f max and of the wavenumber k max at r = r ms are realized 
by the rotation of the geometry. 

The maximum amplitude of the proper growth rate |f max (r)| increases in the similitude 
of the Newtonian angular velocity of a circular orbit £In = 3/4(Mr~ 3 ) 1 / 2 for r — > r ms . 
We show in Fig. 1 the amplitudes of growth rates |f maa ,(r)| for a = M, 0. The ratio 
|f ma x(r)|/(3/4(Mr" 3 ) 1 / 2 ) is restricted to a finite range as pointed by Gammie (2004). Here 
its ratio varies between 1 and 5/(2\/3) « 1.44 when r > r ms and va < 0.1c. The rapidly 
rotating black hole with a = M produces the rapid growth rate which is about one-order 
of magnitude larger than that for a = 0, T max {r- ms ,a = M) = 12.5T max (r ms , a = 0). If the 
marginally bound orbit r mb or photon orbit r ph is realized for a = 0, then the proper growth 
rate |f maa; | reaches to oo. 

It is useful for the magnetohydrodynamics to evaluate the growth rate in MRI in com- 
parison with the speeds of other dynamical processes. The angular velocity of a circular orbit 
Cl is a typical measure to make clear the characteristic speeds of the growth rates, T = T/Cl, 
and to make clear the typical wavelenghth of perturvations, k max = (k ■ va)/&- We show 
in Fig. 3 the dispersion relations with use of f and k for a = M, 0. In the rapidly rotating 
spacetime with a = M, the relative growth rate f is observed to be very large in the vicinity 
of r ms . Its relative value reaches to \T(r r ms )\ ~ 3.75 ~ 5|f (r = oo)|. The unstable range 
of wavenumber k crit is measured to expand in proportion to the shear variable k crit oc % 
On the other hand the relative growth rates f (r) and the critical wavenumbers k crit in the 
static spacetime a = little change at the space r > r ms . These characteristic asymmptotic 
behaiviors of the relative quantities, T max (r — > r ms ) and k cr it, are braught about partially by 
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the variations of Cl(r — > r ms ). When a = and r — > r ms , the circular velocity i>n(r) and the 
proper angular velocity Cl(r) increase in inverse proportion to (r — r^) n , while for a = M, the 
circular velocity v< v \(r) and the proper angular velocity Cl(r) become constant at r ~ 2r ms . 
The rapid rotation of the geometry restricts the increases of vup\(r) and Cl(r) for r — > r ms to 
be V(<p) — vq — = oo — oo ^ 1/2 and f2 = r yQe~ u — •yue~ u = oo — oo — > l/(2\/3). The 
observers in LNRF(ZAMOs) evaluate the growth rate T to be large in comparison with the 
observed angular velocity of a circular motion Q. 

Let's consider the dispersion relation over the extended space around a black hole. The 
rotating black hole distorts the geometry of spacetime not only in the radial direction but also 
in the azimuthal direction. The angular velocity of a circular orbit at a finite height z from 
the equatorial plane, Q(r,z), was represented by Yokosawa, Uematsu & Abe (2004). They 
assumed the pressure force in the vertical direction to be balanced with the gravity. By using 
its angular velocity, Q(r,z), the dispersion relation on the plane, z = contant, is obtained. 
Here we introduce the cylindrical spatial coordinates (w, z, tp) where w and z are defined by 
w — r sin 0,z = r cos 9. The derivative in w or z direction is given by the combination of the 
components of the derivatives in r and 9 directions, e.g., V ro f2 = sin^VJl + cos#V 9 fl The 
distributions of the maximum amplitude of the growth rate, Y 2 max [w ^ z), and the wavenumber 
at the maximum growth rate, k max (zu,z), are depicted in Fig.4, where the Brunt- VdisdYd 
frequency is assumed to be constant, N? = Nj/tl 2 = 0.8 and iV| = O.OliVf, and the 
directions of magnetic field and wavenumber are fixed to be B = B z n^ and k = A^n^. When 
a = M, the amplitude of the maximum growth rate, T 2 nax {w, z), is large over the wide area 
around a black hole. The wavenumber k max {w, z) greatly increases near to the event horizon. 

3.3. Magnetorotational Instability exerted by the Magnetic Fields Induced by 

the Differential Rotation of the Orbiting Frame 

When (5 < 1 and ^vqv^/c 2 ps 1, the term containing v 2 A ^ 2 v^ ip )V^i in the dispersion 
relation (26) exerts on the MRI. When k? = 0, the insertion of this term leads to no change 
of wave mode in the characteristics of perturbations. The sign of this term is contradictory 
to the sign of the change in fluid iV 2 , which means that this term enhances the MRI. The 
dispersion relation for a = M, r = 1.001r ms , z = and kf — is shown in Fig. 5. When 
Alfven velocity va is larger than ~ 0.1c, the growth rate is remarkably enhanced by this 
term and the perturbations with k = become unstable, T 2 (k = 0) < 0. 

When the wave is propagating in the radial direction, 6 oc e r fl*+ i ( r V*+ fc fZ r +fc<p e ) ) waV es 
of perturbations become over stable. The unstable branch of the dispersion relation in phase 
space (kg, k?) is shown in Fig. 7, where the parameters are set to be va = 0.5c, a = M, r = 
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1.001r ms and z — 0. Here the direction of magnetic field is fixed to be B = B e n g . The 
imaginary component of the growth rate IV is also shown in Fig. 7. The wave frequency 
Tj increases in accordance with the radial wavenumber k?. When the Alfven velocity is 
nonrelativistic, va <C c, the unstable branch in the wavenumber space (kg, k?) is restricted 
in a area fcj + kj < k 2 rit (see the left pannel in Fig. 6). Here IV ~ (see the right panel in 
Fig. 6), and then the perturbations exponentially grow. When va > 0.1c, the unstable area 
of the phase space (k e , k?) enlarges to the larger radial wave number, k? > k crit (Fig. 7). The 
oscillating waves exponentially increase their amplitudes over the wide space of wavenumbers 
(k§, kf). 

The above instability is caused by the magnetic induction due to the differential rotation 
of the orbiting frame, 

<9B 

— oc -7 2 ^)M V x B) x V In (Q/a). (34) 
The perturbation of a radial velocity Sv r exponentially increases the magnetic field, 
(91 n B® 

oc -~f 2 v {tp) v n 5v r V r In (SI /a) > 0, when Sv r > 0. (35) 

The increased magnetic field more enhances the motion of fluid according to the equations 
(18)-(20). 



4. Discussion 

Balbus & Hawley (1992) had shown that the maximum growth rate of MRI agrees with 
the local Oort A-value of the disk. It is conjectured that the Oort A-value of a disk is an 
upper bound to the growth rate of any instability feeding upon the free energy of differential 
rotation. They have investigated the reasons behind this remarkably general behavior of 
MRI by using a form of the dynamical Hill equations(Hill 1878). The Hill equations have 
found their widest applicability in the study of collisionless system. We show here that the 
maximum growth rate of MRI in curved spacetime is somewhat different from the general 
relativistic version analogous to the local Oort A-value. 

The local Oort A-value is the growth rate of the dynamical shear instability of a particle 
motion in the field of central potential. Its instability is caused by the displacement of the 
particle orbit in which the angular velocity of a circular orbit keeps a original value, Q = Q . 
When the displacement is positive, i.e. 5r > 0, the centrifugal force is then superior to the 
gravitational force, the radius of the perturbed orbit becomes more larger, and vice versa. 

When the motion is restricted in the equatorial plane, the geodesic equations in LNRF 
are represented with the expressions of the time coordinate, dx^ = df, and the covariant 



- 12 - 



derivative, V^, as 

+ V (r) lna + ^— — V (r) lncj- ( — —J V (r) lne^ = 0, (36) 

j2 (<p) j (p) j (r) _ 

ax + ax ax v (r) lne ^ = 0j (37) 
af 2 af af 

Let's consider the situation of perturbations that the particle is allowed to make small 
excursions from a circular orbit, x^ = e Ml i?o, dx^ = e </3 (f2o ~~ u))dt, by introducing small 
quasi- Cartesian x and y variables: 



W = eP 1 R + x, dx^ =e ip (tt -uj)dt + dy. (38) 



x 



We assume the local force drawing the displaced particle back to its unperturbed lo- 
cation, —Kx, —Ky, where K is the some positive constant. The linearized equations are 
written by 

x — 2Q,?y — 5d x — —Kx, y + Kfi: = —Ky. (39) 

where 

Q- = v H V {r) In - -u w V (r) In x 2 = -2u n ^V (r) In Q and Kf = 2v iv) V (r) In (40) 

and the dots denote the time derivatives, x = d 2 x/df 2 . 

The x and y displacements have solutions to the above equations varying as e* rr , with 

f 4 - f 2 {k 2 - x 2 ) - Kit 2 = 0, (41) 

where f 2 = T 2 — K and k 2 = 2Q?k?. The above dispersion relation leads the following 
expression of the maximum growth rate, 

f 2 = — (AO) 

^ max • V^^J 

Far away from a black hole the functions, x 2 and k 2 , are approximately represented by 

2 dfl 2 2 2ftd(r 2 Q) 

X 2 ~ -— - , k ^— -i. 43 

a In r r ar 

Then the maximum growth rate approaches to the local Oort A-value: 

l r ~l ~ "aSE? (44) 
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Thus the dispersion relation (41) is the general relativistic version analogous to the dynamical 
shear instability in which the maximum growth rate is the local Oort A-value. 

We show the variations of the maximum growth rate r^ ax (r) and the related functions 
R 2 and x 2 for a = M, in Fig. 7. Here we represent the dimensionless variables defined by 

f = f m x = xM « = k/Ci. 

It is a distinctive feature in DSI that at the inner limit, r — > r ms , the maximum growth 
rate |r max (r)| for a = M reaches to oo. Here the function corresponding to the epicyclic 
frequency k approaches to infinitesimal, R — > V2(r — r ms ) — > 0, when r — > r ms while the 
function representing the shear x 2 is finite at the inner limit, x 2 — > 3. The growth rate thus 
becomes |r maa ,| = 3/(2\/2(r/r ms — 1)) — > oo for r — * r ms . The maximum growth rate in DSI 
is remarkably different from that in MRI when a = M and r = r ms . 

This distinctive difference in f max between MRI and DSI is understood with the dis- 
persion relations, (41) and (26). The shear term containing the differential angular velocity 
Vlnf2, i.e., x 2 ', is inserted in the parenthesis at the second term on the left hand of (41) 
while it is in the third term in (26). It leads the different expressions in the maximum growth 
rate between (32) and (42). for r — > r ms , the maximum growth rate T max expressed by (32) 
becomes finite when R 2 — > 0. 

In DSI the term x 2 is produced in the perturbation of "gravitational force" in the 
radial direction (36). In MRI the term containg Vlnf2 appears in the dynamics in the 
azimuthal direction (20) through the perturbation of the azimuthal component of magnetic 
field 8B^. In the Newtonian limit, the angular velocity with the power law dependence 
of radius, f2 oc r~ 3//2 , reduces the differences of the coefficients in the dispersion relation 
between (26) and (41) to disappear, as follows ; k 2 — > fi 2 , k 2 + x 2 4fi 2 , in (26) and 
— X 2 /t 2 — > 4f2 2 in (41). However in rotating spacetime geometry the differential 

rotation of the geometry, t^Vlncj, becomes a main term in these functions. It is concluded 
that there exist remarkable differences in the dispersion relation between MRI and DSI. 

The behavior of maximum growth rate T max (r) in DSI is fairly different from that given 
by Gammie (2004). He showed that the maximum growth rate normalized by Newtonian 
Oort A-value, |r max |/[(3/4)(Mr -3 ) -1 / 2 ], is restricted between 1 and 4/3 for all radii and for 
all a. His method deriving the dispersion relation is different in some aspects from ours. 
He had analyzed the dynamics of perturbations in the Boyer-Lindquist coordinate frames 
(BLCF) in which the observer measures the motion of a particle with clock and rule at 
the distance far from a hole. We have also re-derived the dispersion relation of DSI in 
BLCR with the similar forms as in LNRF, which is shown in Appendix B. Then the function 
corresponding to the epicyclic frequency R reaches to oo when r — > r ms and a — > 1. Thus 
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the maximum growth rate at r — > r ms becomes infinitesimal, T max — > 0, which is contrary 
to the above result in LNRF. 

When a = M, the asymptotic behavior of k? at r — > r ms is extremely different between 
BLCF and LNRF. This difference is caused by the expression of the circular velocity mea- 
sured in each of frame. In BLCF the coefficient Kf in (B9) is described with the velocities, 
vq and v u which become infinity at r — > r ms , while in LNRF Kf in (40) is expressed with v^) 
which is limited within 1/2 for r < r ms . The function Kf in (40) (LNRF) reduces to zero 
for r — > r ms since the covariant derivative lne^ becomes infinitesimal at r — * r ms . The 
dynamics in the azimuthal direction is exerted by the "Coriolis force" KfX. In BLCF the 
rapidly rotating particle perturbed by x is extremely accelerated in the azimuthal direction. 
The dynamics of perturbed motion at r ~ r ms is very different between BLCF and LNRF. 

The perturbed quantities in LNRF are not same as in BLCF. The perturbed velocity in 
the radial direction in LNRF, x, is divided into two components, radial and azimuthal ones, 
in BLCF. Then its azimuthal component becomes not always a first order of small quantity 
in BLCF, and vice versa. Thus the dispersion relation becomes definitely different between 
LNRF and BLCF. 

The LNRF is suitable for the analysis of perturbations around a rotating black hole 
since the dragging of inertial frames becomes so severe in BLCF that the t coordinate basis 
vector (d/dt) goes spacelike at the static limit r . Thus it is concluded that the both types 
of instabilities, DSI and MRI, are remarkably enhanced in the vicinity of a rapidly rotating 
black hole, though their growth rates are not same. 

A summary of the results is as follows. The dispersion relation in MRI is represented 
in general relativistic forms. The growth rate at r ps r ms is realized for a = M by the 
rotation of spacetime geometry. The maximum growth rate |r ma:r | at r = r ms for a = M 
is about one-order of magnitude larger than that for a = 0. When the magnetic field is so 
intensive that va > 0.1c, the over stability mode becomes remarkable and its growth rate 
increases according to the strength of magnetic field. The shear instability, MRI or DSI, is 
largely enhanced by the rapidly rotating black hole though the value of maximum growth 
rate in DSI is fairly different from that in MRI. The dispersion relation at r ~ r ms is suitably 
evaluated in LNRF. In the static geometry the growth rates in MRI and DSI remain almost 
in similar variation as in Newtonian case. 

We are grateful to Professors T. Ishizuka and T. Yoshida for useful discussion and advice 
on the magnetorotational instability in curved spacetime. 
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A. Coefficients of the Affine Connection in the Orbiting Frame 

The rotation one-forms cup are derived from the one-forms uj x shown in (5) (Punsly 
2001; Misner, Thorn & Wheeler 1970). From the one-form u u = (3^dx x , we compute the 
exterior derivative 

dtf = -^dx" A dx x = -u v a A (Al) 
The coefficients of Affine connection Y 1 ^ A are read off the representation u 1 ^ = A uj x . 
Their explicite expressions are as follows: 

r% = e h V f lnd - ln( 7 e^)), 

1% = \l 2 (v n V" Infi - eV)Vnn( 7 e^ M /d), 



T% = e b V y lnd - v^vnYV" ln( 7 e> v) )), 

T % = \v Ql 2 V S \nn - e b v M V § H ie \ v) /a), 

4r =r^=^7 2 V f lnfi, 

ri s =ri f = -^VinQ, 

I* = V*lne^, r^ = -VMne^, 

1^ = - V f ln( 7 e*) + ^ V f In fi, 
= - V § ln( 7 e*) + W n 7 V In n. 



B. The Dynamical Shear Instability in the Boyer-Lindquist Coordinate Frame 

The equation of motion of a particle is described in the Boyer-Lindquist coordinate 
frame(BLCF) by, 

d 2 x x > dx^ dx v 



^ +r ^*=°- (B1 > 

where T x v are the coefficients of Affine connection in BLCF. When the motion is restricted in 
the equatorial plane, these equations are represented by using the expressions of coordinates, 
x ( r \x^\ and the covariant derivative, in LNRF as 



a 2 dt 2 
d 2 x M 



+ V (r) In a + ^^-V (r) ln(cue^) - ( V (r) lne^ = 0, (B2) 
adt V adt / 

{u n [2V (r) lne^ + ^V (r) lnc 
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2V (r) In- 



:i-v; 



;)V (r) In 



to 



X 



dx ir) 



. J adt 



0. 



(B3) 



Let's consider the situation of perturbations that the particle is allowed to make small 
excursions from a circular orbit, r = Rq, ud = Q t, by introducing small quasi- Cartesian x 
and y variables: 

r = Rq + xe~^ , up = tt t + ye^. (B4) 
The linearized equations are then written by 



x - 2v H yV (r) \ne^ + yvj/ {r) \n(uoe^) 



+ e 2 «^V (r) fi 2 V (r) lne^ - ^-V (r) ln(cue 3 ^) 



x = —Kx, 



y + |^[2V (r) lne^ + t; 2 V (r) lncu 

- [2V (r) ln^ - (1 - vl)V {r) In a;] }i = -Xy, 



(B5) 



(B6) 



We follow the convention of writing dots for time derivatives, e.g., x = d 2 x/(a 2 dt 2 ) = 
d 2 x/df 2 . 

The x and y displacements have solutions to the above equations varying as e tTT , with 

-f 2 (k 2 -x 2 + 2K) +K(K -x 2 ) = 0, (B7) 

where 

(B8) 



^ = _ e 2(^)yMr2 2 V (r) lne^-^-V (r) ln(^) 



ft 



V (r) \ne v + -v 2 V {r) In a; 
2 



V (r) ln^-i(l-t; 2 )VM Inc.]} 



x 



V (r) lne^ 



2(n-u) 



V (r) \n(ue^ 



(B9) 



The above dispersion relation gives the maximum growth rate, 

f 2 = ^ 

max — • 



(BIO) 



Far away from a black hole the functions, x 2 an d ft 2 , are approximately represented by 

dtt 2 _„ 2ttd{r 2 Q) 



x 2 ~ 



K 2 ~ 



<i In r r dr 

The above maximum growth rate then has the local Oort A-value: 

1 dVL 



T ~ 

l 1 max | 



2cHnr 



(Bll) 



(B12) 
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We show in Fig. 9 the maximum growth rate f max( r ) an d the related functions, R 2 (r), x 2 (r) 
for a = M. Here we use the angular velocity of a circularly orbiting particle measured in 
LNRF, Cl = (Q — uj)/ a, to define the dimentionless variables such as f = T/Cl, x — x/& an d 
R = R/Cl. The absolute value of the growth rate, |r mox |, for a = M has an upper bound, 
|r maa; | < 1-5, which is a contrast to the maximum growth rate in LNRF. When a = 0, there 
is no large variation in f^ ax (r) for the space, r > r ms . 

At the inner limit, r — > r ms , the corresponding dimensionless functions for a = M 
become 

X 2 -> 3 • 2 8 , R 2 -> cx). (B13) 

The inifinity of k 2 is brought about by i> w — > oo when r — > r ms . However the corresponding 
function in the case of LNRF reduces to infinitesimal at r — > r ms . While the growth rate 
expressed by (BIO) reduces to infinitesimal at r — > r ms , the growth rate expressed by (42) 
reaches to oo. In the curved spacetime the dynamical shear instability in BLCF has the form 
unlike as in LNRF. 
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Fig. 1. — The distributions of the maximum proper growth rate |f max (r)|, proper angular 
velocity Cl(r) and the function 3/4r~ 3 / 2 on the distance from the event horizon. The left 
figure is shown in the case of the extreme Kerr metric and the right one is that of Schwarzscild 
metric. Other parameters are k f = 0, Nj = 0.8fi 2 , N~ = O.OliVf . 
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Fig. 2. — The distributions of the maximum relative growth rate L max {r) an d t ne related 
terms k max (r),x 2 (r),Cl?(r) and Cl shear (r) on the distance from the event horizon: The left 
figure shows the case of the extreme Kerr metric. The right figure shows the case of the 
Schwarzscild metric. Other parameters are same as in Fig.l. 



Fig. 3. — Dispersion relation: The Left figure shows the relative growth rates for the metric 
of an extreme Kerr. The squares of growth rate relative to the angular velocity T 2 = (T/Cl) 2 . 
are depicted. The wavenumber is normalized as k = k • v A /Q. The curves are correspond 
in order to r = 10 2 , 2, 1.5, 1.1, 1.01, 1.001(r ms ) from the top. The right figure shows the 
relative growth rates for the Schwarzscild metric. The curves are correspondent to r = 
10 2 r ms , 1.23r ms , lr ms , 1.2r mb , 1.0r mb , l.2r ph . Other parameters are same as in Fig. 1. 




Fig. 4. — The left figure shows the distribution of the maximum relative growth rate 
r 2 nax (r, for a = M. The right figure shows the distribution of wave number k ma x( r , z) 
at the maximum growth rate. Other parameters are same as in Fig. 1. 
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Fig. 5. — Dispersion relations with the term of magnetic induction due to the differential 
rotation of the orbiting frame. The proper growth rates f 2 are shown as the functions of the 
normalized wavenumber k — k • v^/Q. The curves are corresponding to the values of Alfven 
velocity in unit of light velocity. The metric of spacetime is an extreme Kerr. The position 
is at r = 1.001r ms . Other parameters are same as in Fig. 1. 



- 22 - 



A 




Fig. 6. — Three-dimensional plot of the unstable branch of the dispersion relation with 
va = 1CT 3 . Left figure: The real component of the growth rate T rea i(kg, h?) is plotted in the 
kg(— kgVA/&), K(— kfVA/&) plane. The direction of magnetic field is fixed to be B = B e n§. 
The spacetime of the extreme Kerr is selected and the position is set to be r = 1.001r ms . 
Other parameters are same in Fig. 1. Right figure: Two-dimensional plot of the imaginary 
component of the growth rate Tj m (kg, k?). 




Fig. 7. — The left panel shows the three-dimensional plot of the unstable branch of the 
dispersion relation with v=0.5c. The right panel shows the two-dimensional plot of the 
imaginary component of the growth rate T Im (kg, kf). Others are same in Fig. 6. 
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Fig. 8. — The distributions of the maximum growth rate T max (r) and the related terms 
X 2 (r), R 2 (r) for the dynamical shear instability at the distance from the event horizon: The 
dispersion relation is driven in LNRF. The left figure is shown for a = M and the right one 
is for a = 0. 
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Fig. 9. — The distributions of the maximum growth rate T max (r) and the related func- 
tions X 2 ( r )> ^ 2 ( r )> v uj( r ) m the dynamical shear instability expressed in the Boyer-Lindquist 
coordinate frame. The metric of spacetime is the extreme Kerr. 



